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We discuss the behavior of the Larmor frequency shift and the longitudinal relaxation rate due 
to non-uniform electromagnetic fields on an assembly of spin 1/2 particles, in adiabatic and nona- 
diabatic regimes. We also show some general relations between the various frequency shifts and 
between the frequency shifts and relaxation rates. The remarkable feature of all our results is that 
they are obtained without any specific assumptions on the explicit form of the correlation functions 
of the fields. Hence, we expect that our results are valid for both diffusive and ballistic regimes of 
motion and for arbitrary cell shapes and surface scattering. These results can then be applied to a 
wide variety of realistic systems. 
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I. INTRODUCTION 


The behavior of a system of spins interacting with static and time varying magnetic fields is a very broad topic and 
has been the subject of intense study for decades. A very important application is to the study of spins interacting 
with the randomly fluctuating fields associated with a thermal reservoir. Bloembergen, Purcell and Pound [1] have 
treated this problem using physical arguments based on Fermi’s golden rule and showed that the relaxation induced by 
the fields associated with a thermal reservoir is proportional to the power spectrum of the fluctuating fields evaluated 
at the Larmor frequency, cuo = 7 S 0 , (where 7 is the gyromagnetic ratio and Bq is an applied constant and uniform 
field), which is given by the Fourier transform of the auto-correlation function of these fluctuating fields evaluated 
at Wo- Wangsness and Bloch |5] and then Bloch [3] approached the problem using second order perturbation theory 
applied to the equation of motion of the density matrix and Redfleld, [H (see also m) carried this calculation 
forward to show that the relaxation, indeed, depends on the spectrum of the auto-correlation of the fluctuating fields. 

Another source of randomly fluctuating fields is the stochastic motion of spins (e.g. diffusion) through a region 
with an inhomogeneous magnetic field. To study this problem Torrey [7] introduced a diffusion term into the Bloch 
equation applied to the bulk magnetization of a sample containing many spins (Torrey equation). Cates, Schaeffer 
and Happer |5] then rewrote the Torrey equation to apply to the density matrix and solved this equation to second 
order in the varying fields using an expansion in the eigenfunctions of the diffusion equation. McGregor applied 
the Redfleld theory to this problem using diffusion theory to calculate the auto-correlation function of the fluctuating 
fields experienced by spins diffusing through a (uniform gradient) inhomogeneous field. Recently Golub et al. [TU] 
showed that these two approaches are identical. A useful review of the field is m- 

Another problem that can be treated by these methods is the case of a gas of spins contained in a cell subject to 
inhomogeneous magnetic fields and a strong electric field as in experiments to search for a non-zero electric dipole 
moment (EDM) of neutral particles such as the neutron [12] or various atoms or molecules [13]. This was shown 
by Pendlebury et al. m using a second order perturbation approach to the classical Bloch equation, to lead to an 
unwanted, linear in electric field, frequency shift, (often called a ’false EDM’ effect) which can be the largest systematic 
error in such experiments. Lamoreaux and Golub m showed, using a standard density matrix calculation (Redfleld 
theory), that the ’false EDM’ frequency shift is given, to second order, by certain correlation functions of the fields 
seen by the moving particles. 

Barabanov et. al. m gave analytic expressions for the relevant correlation functions for a gas of particles moving 
in a cylindrical vessel exposed to a magnetic field with a linear gradient along with an electric field. Petukhov et al. 
m and Glayton m showed how to determine the correlation functions for arbitrary geometries and spatial field 
dependence for cases where the diffusion theory applies, while Swank et. al. [19] showed how to calculate the spectra 
of the relevant correlation functions for gases in rectangular cells in magnetic fields of arbitrary position dependence 
even in those cases where the diffusion theory does not apply. Recently Afach et al | 20 ] measured a frequency shift 
that is linearly proportional to an applied electric held (false electric dipole moment) for a system consisting of Hg 
atoms moving in a confined gas exposed to parallel electric and magnetic fields. 

Pignol and Roccia |21| have initiated a program to search for universal expressions giving general results valid 
for all geometries and scattering conditions in the gas and gave such a result for the false EDM effect valid in the 
nonadiabatic (low frequency) limit. Further steps in this direction were taken by Guigue et. al. |22| who provided a 
universal result for frequency shifts induced by inhomogeneous fields in the adiabatic (high frequency) limit and for 
the relaxation rate (Fi) in the case where diffusion theory applies. 

In this work we extend the search for universal expressions of frequency shifts and relaxation for both the adiabatic 
and nonadiabatic (high and low Larmor frequency) limits. 


II. FREQUENCY SHIFTS AND RELAXATION RATES FROM REDFIELD THEORY 


We consider the case of a gas of spin-1/2 particles inside a trap with a gyromagnetic ratio 7 evolving in a slightly 
inhomogeneous magnetic held B{r^ = Bq + b{f). One can define the holding magnetic held Bq = BqBz and the 
Larmor precession frequency ojq = "/Bq. The inhomogeneities b can be taken to have ( 6 ) = 0 where (■ • •) represents 
the ensemble average over all particles in the trap. In addition to this inhomogeneity, the particles can move with a 
velocity v in an electric held E. For simplicity, one can consider that the direction of this electric held is aligned with 
the holding magnetic held: E = Ecz. These particles will experience an effective motional magnetic field E x v/c^. 
The transverse components of the total magnetic inhomogeneity will then depend on the position and the velocity of 


3 


the particles in the trap 


E 

By = by ~;^Vx- (2) 

These transverse inhomogeneities induce a shift 5uj of the precession frequency and a longitudinal relaxation rate 
Ti. Correct to second order in the perturbation, b, the frequency shift Suj, the longitudinal relaxation rate Ti and the 
transverse relaxation rate r 2 , involving the Fourier spectra of the inhomogeneity correlation functions, are given by 
the Redfield theory ElElinil^: 

5uJ — {R,6 [*5'xy(^o) “1“ “1“ *^yy(^o)]} ? 

Fi = 7 ^ {Re + S'yy(wo)] + Im [Syx{uJo) - 'S'a;y(wo)]} 

r2 = y+7"5 ,,(u; = 0 ) 

with 

pOO 

= / e*“"(R,( 0 )R,(r))dr. ( 6 ) 

Jo 

This result is valid in cases where the field fluctuations are stationary in the statistical sense and where the measure¬ 
ments are made over a time scale T 3> Tcon- where the correlation time Tcorr is the time scale for which the correlation 
functions go to zero. In the case of particles evolving in both an inhomogeneous magnetic field and an electric field, 
the frequency shift 5lo and relaxation rate F i can be decomposed as 


(3) 

(4) 

(5) 


5uj = (5wb2 5uje'^ -\- 5ujbe, (7) 

= ri(B2) + + ri(BE), ( 8 ) 


with 


Sujb^ = ^ I™ / e"^‘^'^(bx(0)bx(T) + by(0)by(T)}dT, 

Jo 

poo 

{vx{0)vx{t) + Vy{0)vy{T))dT, 

Jo 


Sujbe =^Re / e^'^^''{bx{0)v^{r)-\-by{0)vy{T))dT, 

Jo 

poo 

ri(B2) =7^ Re/ e'''^°'^{bx{0)b:^{T)-\-by{0)by(r))dT, 

Jo 

poo 

ri(£;2) = ^^^Re / e“°'^(ua;(0)ua;(r) -I-Uy(0)uy(T))dT, 

Jo 

2 

7'i{be) = {bx{0)vx{T) + by{0)vy{T))dT. 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


These Larmor frequency shifts and relaxation rates cannot be further simplified to a form valid for all values of holding 
magnetic field and independent of the particle motion in the trap. However, due to the properties of the Fourier 
transform, there are universal relations that hold for all types of particle motion and all shapes of trap geometry for 
values of Larmor frequency (magnetic field) large and small relative to the inverse transit time of particles across the 
cell, X/v (see below). 


III. SPIN DYNAMICS AND PARTICLE MOTION REGIMES 

In general two length scales describe the motion of a gas of particles in a cell: (i) the mean free path between 
particle collisions noted (ii) the mean distance between two points on the wall which can be evaluated by 

the Clausius expression A = 4VjS where V and S are the volume and the surface of the cell. We define Knudsen’s 
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Figure 1. Classification of the different regimes for a = 5 cm radius spherical cell filled with polarized ®He gas as a function 
of pressure and holding magnetic field. 


number as Kn = y. At high pressure, Kn 1: this is the diffusive regime where the propagation of the particles 
is described by the diffusion equation, characterized by the diffusion coefficient D. At low pressure, Kn 3> 1: this is 
the ballistic regime where the particles travel in straight lines across the cell in free molecular flow. 

The correlation time Tcorr corresponds to the typical time necessary for a particle to probe the magnetic inhomo¬ 
geneity. Since one usually has to deal with large scale inhomogeneities, Tcoti is of the order of the average time between 
successive collisions with the trap walls. Therefore, it depends on the geometry of the trap and on the properties of 
the particle motion inside this trap. In the case of a gas at atmospheric pressure, the correlation time is about 1 s 
for a cubic trap with 10 cm sides. For rarefied gas confined in the trap, Tcon- is approximately equal to 1 ms. This 
time scale can be compared with the Larmor precession frequency wq- The limit when ojq is much bigger than 1/Tcorr, 
is called adiabatic regime. This regime can be interpreted as the particle spins following the local magnetic field. It 
is also valid when the particles are moving slowly in the trap or if they encounter a great number of collisions with 
other particles between two collisions with walls. In contrast, the regime is called nonadiabatic if loqTc ^ 1. This 
limit physically appears when the particles are able to probe the whole magnetic inhomogeneity within times shorter 
than a Larmor period. It is also sometimes refered to as the regime of motional narrowing. This phenomenon can be 
observed in systems immersed in very weak magnetic fields or if the thermal particles are in a ballistic regime in a 
small container. 


For a given trap geometry, these regimes depend on the pressure of the spin gas and on the holding magnetic field. 
Fig. 0 shows this classification as a function of pressure and holding field for a ^He gas contained in a spherical cell 
with 5 cm radius. The super-adiabatic regime corresponds to the situation where the gas is in a diffusive regime and 
the spin motion is adiabatic between two interparticles collisions. In this case, we have the condition WoTcqu ^ 1, 
with Tcoii the time between two interparticles collisions. The correlation functions calculated in m are valid in this 
region as is Eq. (12) in |3]. 


To illustrate this classification, let us consider some realistic systems of particular relevance for our study. The 
cylindrical RAL/Sussex/ILL trap El, used to measure the neutron EDM, contains ultracold neutrons (UCN) and 
a mercury comagnetometer, immersed in a 1 fj,T holding magnetic field. The small particles number of each species 
and the size of the trap (47 cm diameter and 12 cm height) lead to a large Knudsen’s number and so to the ballistic 
regime. However, the speeds of the two particle species are very different: while the mercury atoms are at thermal 
equilibrium and have an average speed of several hundred meters per second, the UCN are moving at a few meters 
per second. Therefore the mercury comagnetometer is in the ballistic nonadiabatic regime whereas UCN are in the 
ballistic adiabatic limit. In the case of a gas at atmospheric pressure, such as a polarized ^He gas HU, in a several fiT 
holding field, the particle motion follows the diffusion equation and the number of Larmor precessions done by the 
spins between two collisions with the walls is very high. This kind of systems is thus in diffusive adiabatic regime. 


As shown in ED 1221, the leading order of frequency shifts (|^ , and ( [TI| ) for adiabatic and nonadiabatic regimes 
can be expressed as powers of ujQTcori or l/woTcorr- To do so, we apply a succession of integrations by parts of the 
integrals defining the frequency shifts. This is the purpose of the two next sections. The first one presents the 
simplified expression of the frequency shift in the adiabatic regime. The nonadiabatic regime will be considered in 
the second next section. The well-known case of uniform magnetic gradients will be discussed in Section VI. 
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IV. ADIABATIC REGIME: HIGH MAGNETIC FIELD OR SLOW PARTICLES, ARBITRARY FIELDS 


The adiabatic regime corresponds to systems which satisfy uiq ^ 1/Tcorr- We want to expand the frequency shift 
Q in power series of 1/wo- way to obtain such an expression consists in applying several integrations by parts 


POO 

/ sin(a;oT)/(r)dT = 

Jo 


— cos{ujot) 


Wq 


/(t) 


1 f°° df 

H-/ cos(wot) —(T)dr, 

<^0 JO 


(15) 


cos(a;or)/(r)dT = 


sin(a;oT) 


Wo 


/M 


_ 0 ^0 Jo 


- rsin(wor)^(r)dr. 
Jo Jo dr 


(16) 


Equation ( |15[ ) and (|16| assume that the function / and its derivative are integrable. We denote / = ^. Using the 
fact that t 


re correlation functions go to zero at infinite time, we can write 

2 2 

+^y) - + by iO)bym 


^ J COs{uJoT){ba;{0)bx{T) +by{0)byiT))dT, 


(17) 


_ {Vx + w) 


2 c"‘wo 


2 j^2 


+ 


YE 

2c‘^ujt 


0 Jo 


cos(a;or)(ua;(0)ux(r) + Uy (0)uj,(r))dr, 


(18) 


Sujbe = 


Ye 

Ye 


{bxiO)VxiO) + byiO)VyiO)) 


cY Jo 


cos{u!oT){bxiO)vxiO) + 6 y( 0 )i;y( 0 ))dT. 


(19) 


Making the reasonable assumption that the correlation functions and their derivatives are continuously decaying to 
0 for r —> oo, we apply the Riemann-Lebesgue lemma |23| and arrive at the conclusion that the last terms in Eq. (171, 


(181, (191 go to zero faster than the other terms in each equation. The frequency shift expressions can be written as 


(see also Eq. (27), (28), (291): 


2 2 

< 5 wb 2 = ^{bl + bl) - ^{bx{d)bxY) + by{Q)by{0)) + O {I/YoTcoY^) . 


( 20 ) 


{vl +vl) +0 (l/(woTcorr)^) , 


( 21 ) 


5ujbe = J^i^{bx{0)vx{d) + 6y(0)uj^(0)) + O {l / YqTcotY) ■ (22) 

C Wq 

Using the expressions for the derivatives of the correlation functions presented in the Appendix and assuming that 
velocities in different directions are uncorrelated and (u^) = {Vy) = {vD = we obtain 

2 2 

^ ^ ^ (l/(jJ0Tcorr)') , (23) 


A / 2^ 


■ O (l/((^0'^corr) 


(24) 
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= 


c^ujI 


,dbx 2 \ 


,dh 

'l>y 


+0 (l/(a;orcorr)‘^ 


(25) 


These results are presented in Table The first term in Eq. (23) corresponds to the leading order of the frequency 
shift in adiabatic regime |22| . 

It is instructive to note that these results can be obtained in another way. One can rewrite Eq. ([^ 


nOO poo 

S,,{oj)= E^^B,{0)B,{T))dT = e“V(^)d7 

Jo Jo 

and expand /(t) in a Taylor series 


(26) 


Sij{uj) = 


= /( 0 ) + /( 0 ) 


/(O) + /(0)t + • • ■ + /^”H0)— j 

d /("Ho) 


d{iuj) 




= /(O)- - /( 0)4 + /(O) A + • • • + 


'0 

(- 1 )” 


.n+l^'n—1 


(27) 


Taking the real part (that is, the relaxation rate for the B ^, terms and the frequency shift for the EB term, we 
obtain 


1 


Re [S., (c.)] = -/(O) ^ + ... , (28) 

which is equivalent to Eq. (14) of [55] and ( |T^ above but now we have another form of the correction term and we 
can use this to calculate the frequency range where the first term is a good approximation. For the imaginary part 
{B^,E^ frequency shifts and EB relaxation) we find 


Im [^^(w)] = 


/(O) /(O) 




LO UJ'^ 

which is equivalent to Eq. (17) of the same paper and Eq. (|20| above. In addition we find. 


1(B2) = 7 


1 


1 


-2 (^-(O)^-(O) + by{0)by{0)) + ^(&,.(0) 6,(0) + by{0)bym 


Wn 


Wn 


1 (B 2 ) 


1{BE) 




2'y^E 


1 


1 


- 2 + T(^^x( 0 )n,( 0 ) +ny( 0 )ny( 0 )) 


Wn 


OJn 


1 


1 


— (6,(0)n,(r) + by{0)vy{T))\r=o + ^(6,(0)i;,(0) + by{0)vy{0)) 


Wo 


Wn 


Using the results presented in the Appendix, we can write 


1 (B 2 ) 


2 w 2 


^ V 


{bl +bl)) + O ^1/ (woTcorr)"') , 


— 2 '^{'^x'bx T '^y'by} T O (l/(woTcorr) 

CJnC 


YE 

4 

2YE 


bl(S£:) — 2 {bx'^x T byVy) -f O (l/(wo'7"corr) ) 


WqC 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


and ([Ml as volume averages of the velocity 


We expressed the correlation function derivatives in (231, (241, 
and magnetic field, in the adiabatic limit. These expressions are therefore independent of the particle motion in the 
cell. The high frequency limits for Sugi, and r 2 (B£;) are universal. 

The first term in (331 behaves as I/wq and has been calculated in in the diffusive adiabatic regime. However the 
diffusion theory breaks down at times shorter than the collision time Tcoh, so that at high frequencies, (woTcoh ^ 1) 
the spectrum deviates from that expected on the basis of diffusion theory. The high frequency (super-adiabatic) limit 
is correctly given by uni, using a correlation function that is valid for all times, namely the first term in ( |33| ) goes 
to zero as the velocity is initially uncorrelated with position, and the very high frequency behavior goes as (l/wg). 
The result of m shows how the behavior goes from the (~ 1 /wq) predicted by diffusion theory at high frequencies 
to (~ I/wq) as woTcoii becomes on the order of or greater than 1. 
























V. NONADIABATIC REGIME: WEAK MAGNETIC FIELD OR FAST PARTICLES, ARBITRARY 

FIELDS 
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We now consider the nonadiabatic limit woTcorr ^ 1. To expand the frequency shifts expressions in terms of power 
of (jJo, we simply apply the same procedure of recursive integrations by parts, changing the part that is integrated: 


sin(a;or)/(r)dr = 


sin(a;or) / f(t)dt 


— Wo / cos(a;or) / /(t)dtdr 


(36) 


J 0 


POO r PT 1 OO noo PT 

/ cos(a;oT)/(T)dT = cos(a;o'r) / f(t)dt + wq / sin(woT) / /(t)dMr. 

Jo . Jo . 0 Jo Jo 


When applying these relations to Eq. (§, ( [To| and ( [Tl] ), we obtain: 


Suje^ = — 

+ 


2c* 

2c* 


uoix'^ + y^) 

POO 

Wq / sin(wor)(a;(0)a:(r)+ y(0)y(T))dr 


(37) 


(38) 


StORE = — 




{b^x + byy) 

pOO 

sin(woT)(6a;(0)a:(r) + by{0)y{T))dT. 


c 2 ^ 

f 


(39) 


One can see that the last terms on the right-hand sides of Eq. (381 and (391, as well as Eq. ([^ involve Fourier 
transforms of correlation functions which depends exclusively on position. Since we are in the limit woTcon- ^ 1 we 
need to consider only the first order expansion of the involved trigonometric functions: sin(a;o'r) « ojqt. (For times 


> 


the correlation function goes to zero.) Applying this method to Eq. ( [M| and ( [M| ), we obtain 


\ 

5uJb2 ~ -—UJq 


r{b:^{0)bx{r) + &y( 0 ) 6 y(r))dT 


(40) 


dujE^ « -^^uJo{x +y ) 


72^2 


-Wn 


r(a;(0)a:(T) -f y{0)y{T))dT 


(41) 


Sujbe ~ — 


72 ^ 


{bj:X + byy) 


(42) 


ye 


-{bxiO)x{T) + by{0)y{T))dT. 


The last terms in Eq. (411 and (421 can not be calculated for any arbitrary trap geometry. But one can see that they 
behave as uigTcorr when woTcorr goes to zero. This means that the expressions of the frequency shifts are dominated 
by the first term on the right hand side. These results are presented in Table |T] 

Similarly, since coswqt « 1, Eq. (121, (131 and (|l4|)) become 


ri(B 2 ) =7 / 

Jo 

Ye^ 


ri(£) 2 ) = 


l(BE) 


2YE 


{bx{0)bxY) + by{0)by{T))dT^ 


Y / {x{0)x{t) + y{0)y{T))dT, 

Jo 

POO 

Wo/ {bxiO)x{T) + by{0)y{T))dT, 

Jo 


(43) 

(44) 

(45) 


from which the low frequency limits follow immediately. 

























Frequency 

Adiabatic 


Nonadiabatic 

shift 

(UCNs) 


(Hg) 

Sujgi 

&,{bl+bl} + ^y}{\Vbx 

^+\y?) 

^^0 r{hz{f))bz(T) + by{Q)by{T))AT 




/ 2 I 2\ 

Sujbe 

-y.[yvl) + yvl)) 


y{bzX + byy} 


Table I. Expressions of the leading terms of the frequency shifts induced by the transverse magnetic and motional fields, in the 
adiabatic and nonadiabatic limits. 

VI. MAGNETIC FIELD LINEARLY DEPENDENT ON POSITION (UNIFORM GRADIENTS) 

In this case, which is what is usually treated theoretically and applies to most experimental situations, one can 
simplify the derivatives in terms of trajectory correlation functions without any evolution equation and derive a variety 
of relationships. Let us consider a magnetic inhomogeneity b dependent linearly on the position of a spin: 

K = G^x, (46) 

by = Gyy, (47) 

where the relation Gx + Gy = = —Gz holds by the divergence theorem. 


A. General relations for fields with linear gradients and cylindrical symmetry 

1. Expressions relating frequency shifts with frequency shifts 

In the common case of cylindrical field symmetry, (with arbitrary cell shape) the correlation functions of interest 
are (ua;(0)va;(r) + Uy(0)vy(r)), (a:(0)ua;(r) +y{0)vy{T)), {x{0)x{t) + y{0)y{T)) which satisfy the relations 


(x(0)ua:(T) + y{0)vy{T)) = - — {x{0)x{t) + y{0)y{T)), 

K(0)u„(r) + Uy(0)uy(T)) = ^{x{Q)vz,{t) + y(0)uy(r)) = --^{x{0)x{t) + 2/(0)?/(t)). 


Then 


Sxv (w) = -iujSxx (w) + (a;^ + , 

Svv (w) = -iuiSx^v (w). 


According to Eq. ®[T^[TT]) we find 


5ujbe = KBE^e{Sxv{i^)) = -ArB£:Im(5'„„(a;)/a;) = - 


Kbe Slob^ 
Ke2 uj 


(48) 

(49) 

(50) 

(51) 

(52) 


with Kbe = = s “ ’ ~ 2 c^- expression was obtained in [14] for the case of particles 

moving in a cylinder with specular reflecting walls and no gas collisions, but our result holds for any cell shape and 
type of particle motion. Equation (391 can be written as 

l^EGz 


Sujbr = 


2 c 2 




YEGz 

2 c 2 


Wo / sin(woT)(x(0)x(T) + y(0)y(r))(ir 


j'^EGz, 1 2\ 4ifwo 


2 c2 


-kx^Ey^)- 


- Lt? 


6ojb^, 


(53) 


(54) 


so that Eq. (541 represents a method of measuring the linear in E shift without applying an electric field. To do this 
one would apply a known constant gradient and look for a frequency shift dependent on the square of the gradient. 
The possibility of the volume average held being changed by application of the gradient can be accounted for by taking 
the part of the shift proportional to the square of the gradient. Another method would be to measure the relaxation 
rate due to then application of the gradient as discussed in the next paragraph. 
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2. Expressions relating frequency shifts and relaxation rates 


As the correlation functions defined in Eqs. (|^- [l4|) are all causal, that is, they are zero for r < 0, their real and 
imaginary parts are related by a dispersion relation |24| and we can write 

5ujbe = Kbe [{x^ + y^) - wim (w)]] 

Re [Sxx (w')] 


= Kbe 


= Kbe 




du)' 


= K 


BE 


= K 


BE 


{x^ + y^)- 


{x^ + y^)- 


{x^ + y^)- 


1 


TT 2Kb2 

_J_uP 

Kb^ tt 
1 u;^ 

Ke^ tt 


UJ — UJ 

pOO Y 


(55) 

(56) 


' —<50 
pOO 


' —<50 
/*<50 


LU — LU' 

oP- — 

ri(_E2) (uj') 
/_oo (w^ -uj'^)uj' 


duj' 


fduj' 


(57) 

(58) 

(59) 


Equations (58 1 and (59) are particularly interesting because they allow measurement of the frequency dependence 
of Sujbe, the shift, linear in E, that produces a serious systematic error in the searches for particle electric dipole 
moments without application of an electric field. By applying a gradient, larger than any existing gradients 
and measuring ri(B 2 ) (w) one can reconstruct the frequency dependence of Sujbe- For the case of a non-cylindrically 
symmetric cell we can apply relatively large gradients dbx,y/dx,y and thus measure, separately, the spectra of the 
correlation functions in the two directions. While according to ( |58[ ) we need to know the relaxation for all frequencies, 
the necessary range of measurement is limited because the known high and low frequency limits are reached rather 
quickly ( 30|43 1 . Substituting (13 1 into (59 1 we obtain a form of the relation that has been obtained by another method 


m- 


3. Expressions relating relaxation rates with relaxation rates. 

For completeness we give relations between the relaxation rates which are abtained in a similar way: 


1 (B 2 ) 


Ke- 

Kb2 




1 (B 2 ) 


Ke^ 

Kbe 


ujV 


1{BE)- 


(60) 


The relaxation caused by the electric field alone, ri(£; 2 p has been discussed in 


B. Adiabatic regime: high magnetic field or slow particles for fields with uniform gradients 


We specify Eq. (23), (24), (25) to the case of a uniform gradient: 

^XlB^ = + ^y) + EES + Gl{Vy)} + 0(l/a;oT®orr) 


0 


2loq 


^W_E2 = + vl) + 


Sujbe = + 0{1/ujIt^„„). 

C UJr. 


Similarly, the relaxation rates can be expressed as: 

1 


Fqs^) = -7^^ [Glixvx) + Gliyvy)] + 0(l/a;^r^^„J, 

UJq 


r - 

^ i{BE) - 


1 

Wo 


{GxXVx + Gyyvy) + 0(l/a;gr^V,) 


(61) 

(62) 

(63) 

(64) 

(65) 


The first term in Eq. (611 corresponds to Eq. (18) in |22]. It is remarkable that it is possible to derive a simple and 
universal expression for the third order term cx in Eq. (611. To our knowledge this third order term has never 
been calculated before. 
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C. Low field, high velocity limit for fields with uniform gradients 


For uniform gradients (471, the expression of Sujbe (421 can be simplified to 


Sujbe = + Gyy'^) + [ T{Ga:x{0)x{T) + Gyy{0)y{T))d7 

c C Jo 


-{GxX d-Gyy ) + 0 {ujqT^^^^). 


( 66 ) 

(67) 


Also (40) becomes 


-7^ 

(5a;s2 « y Wo r(G^a;(0)a;(r) + Gy?/(0)7/(T))dT 


( 68 ) 


and (41) 


Slue^ = — 


72^2 

2c4 


/*oo 

uJo{x^ + 2 /^) + t{x{0)x{t) + y{0)y{T))dT. 


(69) 


Similarly 


pOO 

ri(B 2 ) = 7 M {GIx{0)x{t) + Gly{0)y{T))dT, 
Jo 


p_ , 


UE^)- ,4 


r 


i(B£;) 


2 p2 roo 

Lul / (a:(0)a;(r) + y(0)y(r))dT, 

Jo 

2j-^E 


n PCX) 

-Wo / (Ga;a;(0)x(r) -f Gyy(0)i/(r))dT, 

Jo 


(70) 

(71) 

(72) 


so that in the low frequency limit there are no universal or quasi-universal expressions. 


VII. CONCLUSION 


In this paper we have investigated the asymptotic behavior of the spin-relaxation and related frequency shifts due to 
the restricted motion of particles in non-uniform magnetic and electric fields. Simple universal expressions (valid for 
any form of gas container and any spatial form of the field) were obtained for the observables 6uj and F 1 for adiabatic 
and nonadiabatic regimes of spin - motion. The remarkable feature of all our results is that they were obtained without 
any specific assumptions about the explicit form of the correlation functions. Hence, we expect that our results are 
valid for both diffusive and ballistic regimes of motion. These results can then be applied to a wide variety of realistic 
systems. They are especially important in the context of experiments searching for the electric dipole moment using 
trapped particles, for the frequency shifts proportional to electric fields are of utmost importance. In particular we 
have given general relations between various frequency shifts and relaxation rates. 
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Appendix A: Useful expressions of correlation functions 


We derive three useful relations for the derivative of certain correlation functions in terms of volume averages. Let 
us consider an inhomogeneity bi with i = x or y. 


^(6»(0)&i(r))|^=o = {bi{0)b,{0)) 


= {h^V,) + {b.^Vy) + {b,—V.) 

= l{v-yb^) 


(Al) 

(A2) 

(A3) 

(A4) 


dr^ 


{h{0MT))\r=O = {briO)hm 

= -ihioMo)) 

.-((vvLT) 


dbj 

dx 


dbi 

dy 


= /fA 


dbi 


dbj 

dz 




dbi 


(A5) 

(A6) 

(A7) 

(AS) 

(A9) 


— (fo*(0)nj(r))|^=o = -(&i(0)vj(0)) = -{{Vh ■ v)vj) 

/dbi 2\ 


(AlO) 

(All) 
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